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ABSTRACT 

Emission line fluxes from cool stars are widely used to establish an apparent emis- 
sion measure distribution, EmdupplTc), between temperatures characteristic of the 
low transition region and the low corona. The true emission measure distribution, 
Emdt{Te), is determined by the energy balance and geometry adopted and, with a 
numerical model, can be used to predict E md (Tc) , to guide further modelling. The 
scaling laws that exist between coronal parameters arise from the dimensions of the 
terms in the energy balance equation. Here, analytical approximations to numerical 
solutions for Emdt{Tc) are presented, which show how the constants in the coronal 
scaling laws are determined. The apparent emission measure distributions show a min- 
imum value at some To and a maximum at the mean coronal temperature Tc (although 
in some stars, emission from active regions can contribute). It is shown that, for the en- 
ergy balance and geometry adopted, the analytical values of the emission measure and 
electron pressure at To and Tc, depend on only three parameters: the stellar surface 
gravity and the values of To and Tc. The results are tested against full numerical solu- 
tions for e Eri (K2 V) and are applied to Procyon {a CMi; F5 IV/V). The analytical 
approximations can be used to restrict the required range of full numerical solutions, 
to check the assumed geometry and to show where the adopted energy balance may 
not be appropriate. 
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1 INTRODUCTION 

The observed fluxes of emission lines from cool stars, in- 
cluding the Sun, are now routinely used to derive appar- 
ent emission measures (-Bmapp) for given lines and apparent 
emission measure distributions [i?mdapp(rc)], by using dif- 
ferent lines formed over a range of electron temperatures 
(Te). The precise definition of the emission measure differs 
between authors and that adopted here is given in Section 
2. 

When the International Ultraviolet Explorer (lUE) was 
operating, the _Emdapp(re) of cool stars was known only be- 
tween Tc ~ 10"' K and ~ 2 x 10^ K (see e.g. Jordan et 
al. 1987). Early X-ray satellites and the Extreme Ultravio- 
let Explorer (EUVE) together provided information above 
~ 8 X 10^ K, and under favourable circumstances the EUVE 
could detect a few lines formed around 3 x 10^ K (see e.g. 
Drake, Laming & Widing 1995). The Goddard High Resolu- 
tion Spectrograph (GHRS) and the Space Telescope Imaging 



Spectrograph (STIS) on the Hubble Space Telescope (HST) 
have given improved spectral resolution and sensitivity in 
the ultraviolet regions (e.g. Wood et al. 1996). In particular, 
in addition to lines formed in the low to mid transition re- 
gion (at 10* - 2 X 10^ K), the spectrum of an M-dwarf flare, 
obtained with the GHRS, showed the forbidden line of Fe xxi 
at 1354 A (Maran et al. 1994). Also, spectra of G/K dwarfs 
obtained with the STIS showed the forbidden lines of Fe xil 
at 1242 and 1349 A (Jordan et al. 2001), formed in their 
upper transition region or inner corona. The Far Ultraviolet 
Spectroscopic Explorer (FUSE) has observed the resonance 
lines of O vi (formed around 3 x 10^ K), a number of lines 
that are formed below 10^ K and further forbidden lines of 
Fexviii and Fexix (Redfield et al. 2003). Both Chandra and 
XMM-Newton now provide extensive information from lines 
formed at Tb > 8 x 10^ K (see e.g. Sanz-Forcada, Favata 
& Micela 2004). Thus the Emdi,pp{Tc) of a range of stars 
is now reasonably well constrained over all the temperature 
range of interest, and in some stars the observations of for- 
bidden lines of highly ionized iron with the STIS and FUSE 
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provide simultaneous measurements in the lower and upper 
transition region/corona. 

The true emission measure distribution [Emdt{Te)] is 
determined by the energy balance in the outer atmosphere 
and the actual geometry. By making numerical calculations 
in a chosen geometry, theoretical values of both Emdt{Te) 
and Emda,pp{Te) can be predicted and the latter can then 
be compared with that observed. In an example of this 
approach, Sim & Jordan (2003a) used EUVE observations 
of e Eri (HD22049, K2 V) to determine EnuLpp{Tc). The 
numerical model assumed that in the upper transition re- 
gion there is an energy balance between the divergence 
of the conductive flux and the radiation losses and, as a 
boundary condition, was required to match the observed 
minimum in i5mrfapp(To). A spherically symmetric geom- 
etry was adopted and a non-thermal pressure term was in- 
cluded in the equation of hydrostatic equilibrium. The ob- 
served EmdiippiTc) (and observed electron densities, Nc) 
could be fitted by the values from the model provided the 
area filling factor was ~ 1 in the inner corona, and ~ 0.2 
in the mid transition region. More recently, Ness & Jor- 
dan (2008) have included an analysis of the X-ray spec- 
trum obtained with the Low Energy Transmission Grating 
Spectrometer (LETGS) on Chandra to obtain an improved 
observed Emdapp{Te) and also the relative element abun- 
dances. The atomic data from CHIANTI (v 5.2) (Landi et 
al. 2006) were used to update the values of the electron pres- 
sure (Pc). Improved area factors, which were slightly smaller, 
were also found by iterating the first solution with a variable 
area factor. 

Early work on the implied area of emitting material in 
the solar transition region was carried out by Kopp (1972). 
Also, Torriccelli-Ciamponi, Einaudi &: Chiuderi (1982) used 
the presence of a minimum in the Emdapp{Te) to constrain 
models for the heating of solar loop structures, in an ex- 
tension of the approach used by Rosner, Tucker & Vaiana 
(1978). 

Hearn (1975, 1977) derived scaling relations for the en- 
ergy losses by conduction and radiation from coronae, by 
assuming minimum energy loss from a corona. His approach 
has some similarities to that adopted below and we compare 
his predictions with our results in Section 2. 

The results of a range of numerical models made for 
£ Eri by Ness & Jordan (2008) showed that, provided the 
heating function can be expressed as an energy flux, there 
exist scaling laws between the calculated coronal emission 
measures, electron pressures, temperatures and the value of 
the stellar gravity, g{rc), at the radial distance at which the 
coronal temperature, Tc, is reached. This is to be expected, 
since these scaling laws can be derived from simple dimen- 
sional arguments, with the constants being determined by 
the specific assumptions made. However, in numerical cal- 
culations in other than plane parallel geometry, and in cal- 
culations that include the non-thermal pressure term, it is 
difficult to see exactly what determines the constants. 

Here we adopt some simplifications regarding the varia- 
tion of Pe and the geometry, to find which chosen parameters 
control the full solutions to the Emdt{Te). It is shown that 
the analytical solutions given below depend only on the stel- 
lar surface gravity (p,) and the boundary temperatures, Tc 
and To. The values of Pp, Pc, Emt{To) and Emt{Tc) are all 
determined by the choice of g*, To and Tc. 



The simple approach adopted is also useful as a starting 
point for testing the following questions: does the emission 
have to come from some fraction of the surface area; is local 
deposition of a heating flux (other than from conduction) 
required below the corona; is a mean corona, rather than 
closed magnetic structures (with heights smaller than the 
isothermal pressure-squared scale-height) appropriate? 

Although To is similar in the main-sequence stars stud- 
ied (~ 2x 10^ K), it is significantly higher in the F-star Pro- 
cyon (HD61421. q: CMi F5 IV/V) and other evolved stars, 
including single giants and RS CVn binaries (as pointed out 
for Capella by Dupree et al. 1993). The giant stars and 
RS CVn binaries arc not considered here because we as- 
sume a plane parallel geometry, which becomes inappropri- 
ate in lower gravity stars and in the coronae of close binaries. 
Also, consideration of their escape velocities suggests that 
the highest temperature material is likely to be magnetically 
confined, rather than occurring in a quiescent corona with 
a large scale height. 

The simple theoretical model adopted for stars with a 
quiescent corona is set out in Section 2, together with the an- 
alytical relations derived. Comparisons between the simple 
analytical results and those derived from the full computa- 
tional models for e Eri, by Ness & Jordan (2008), are made 
in Section 3. Procyon is used as a further example in Section 
4. The conclusions are discussed and summarized in Section 
5. 



2 THEORY: STARS WITH A QUIESCENT 
CORONA 

In the Sun, the maximum in the observed Erndupp^Tc) 
occurs at the average temperature of the inner quiescent 
corona (Tc ~ 1.6 x 10® K). Except in observations at or 
above the limb, the emission lines are formed predominantly 
over the first pressure-squared isothermal scale height, since 
the lines are coUisionally excited, with fiuxcs oc TV^ (Sec Sec- 
tion 2.2.) As apparent in movies made using Yohkoh data, 
the emission from the inner corona, away from active re- 
gions, shows little inhomogeneity and very few variations 
with time. If active regions are present, the emitting mate- 
rial is observed as a gradual decrease in the Ernds.pp(Tc) at 
temperatures larger than Tc. In very early spatially unre- 
solved studies, Neupert (1965) observed the variation with 
time of lines formed at different values of Tc, over more than 
one solar rotation period. These support the above picture. 

In other main-sequence stars we expect a quiescent 
corona to be present, as well as active regions. Tc is as- 
sumed to correspond to the temperature at which the ob- 
served EmduppiTc) has its maximum value. At worst, this 
assumption gives an upper limit to Tc. Ideally, the contribu- 
tion of active regions should be studied through rotational 
modulation of high temperature X-ray lines, but this is cur- 
rently diSicult, owing to the large amount of observing time 
required. 

2.1 The upper transition region 

The full numerical calculations were made using a spher- 
ically symmetric atmosphere and included a non-thermal 
pressure term in the equation of hydrostatic equilibrium. 
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The relevant equations have been given in Sim & Jordan 
(2003a) and/or Griffiths & Jordan (1998). In the analyti- 
cal results below, we adopt a plane parallel geometry and a 
constant emitting area. In Section 3 we compare the analyt- 
ical results with those from the full spherically symmetric 
numerical models calculated by Ness & Jordan (2008). 

The true emission measure for a given line is defined as 



£mt(0.3) = [ Nn 

Jsh 



N^dh 



(1) 



where h is the height, 5h corresponds to the region over 
which logTe changes by ±0.15 about the optimum tempera^ 
turo of formation for the line and Nh is the number density 
of hydrogen. 

In plane parallel geometry, the apparent emission mea- 
sure for an optically thin line is given by 



Smapp(0.3) = i / NuNedh , 



(2) 



where the factor of 1/2 allows for the photons emitted in 

the outwards direction. 

Eqn. (1) can also be expressed as 



Smt(0.3) 



0.86Pe dh 
V2Te dTe 



(3) 



where the electron pressure is defined as Pe = NgTe, and 
the factor of 0.86 arises from A'h — 0.86A'o. Here it has 
been assumed that and dh/dTe are constant over the re- 
gion in which an individual line is formed. The former is a 
good assumption, but the latter can initially be less accu- 
rate. However, in actual calculations of line fluxes using a 
theoretical model, any variation of dh/dTc over the region 
of line formation is taken into account. The starting values 
of i5mapp(0.3) are used to find the initial _Bmdapp(0.3)(Tc), 
but following the calculation of the line fluxes, including the 
full contribution function for each line, iJm(iapp(0.3)(Tc) is 
optimised. Relative clement abundances are also adjusted 
during this process. The details are given in Ness & Jordan 
(2008). Note that eqn. (3) gives an expression for the true 
emission measure for a line formed at a given Tg. Eqn. (3) 
is also used to define the true emission measure distribu- 
tion, Emdt{0.3){Te), since the quantities on the right-hand 
side are all differentiable functions of Tg. For simplicity of 
presentation, the label 0.3 is omitted in the equations below. 

The energy flux carried by thermal conduction, Fc{Te), 
is given by 



Fc(Te) 



5/2 



dTe 

dh 



(4) 



5/2 

where kTc is the coefficient of thermal conduction 
and, from Spitzer (1956), k is taken to be 1.1 x 
10~ erg cm' s"^ K"''/^. Here, as in our full numerical 
calculations, we ignore the small variation in k with A^c and 
To (Spitzer 1956). In e Eri, this amounts to only 30 per cent 
between log To = 5.3 and 6.5. There will be a difference be- 
tween Fc(Te) in the analytical and numerical, spherically 
symmetric, solutions, but given the limited extent of the re- 
gion considered, this is not expected to be large in dwarf 
stars [see comments after eqn. (12)]. 

Eqn. (3) allows Fc{Te) to be expressed in terms of Errit. 

I.e., 



0.86k PfT, 



.3/2 



(5) 



y/2 Emdt(Te) ' 

From eqn. (5) it is clear that one cannot use a boundary con- 
dition that sets Fc(Te) to exactly zero at some base temper- 
ature, since this would lead to an inflnite value of Emdt{Te). 

As carried out in earlier work [see Jordan & Brown 
(1981)], eqn. (5) can be differentiated to give 

d\og EmdtjT,) ^ 2dlogPe _ dlog(-fc(re)) , . 

dlogTe ' dlogTe dlogTe ' 

So far, no assumptions have been made about the en- 
ergy balance. It is now assumed that the corona is heated 
by a flux of energy from lower layers and that this energy 
is not dissipated until high in the corona, far above the first 
pressure-squared isothermal scale-height, H, over which the 
spatially averaged stellar emission lines arc mainly formed. 
(The same situation is relevant to solar lines observed near 
Sun-centre.) In this case, the heated region can be studied 
only through solar observations above the limb. In the Sun, 
Te rises slowly up to a height of 0.70 Rq, far higher than 
H ~ O.O6-R0, while the non-thermal velocities (interpreted 
using Te, rather than the unknown ion temperatures Ti), 
continue to increase (Landi, Feldman & Doschek, 2006). 

Thus, below the heated region, it can be assumed that 
the divergence of the thermal conductive flux is balanced 
by the radiation losses. If am^ dissipation of the mechanical 
heating flux (-Fm) were present in the upper transition region 
and inner corona, this would add a term —dFm/dh to the 
right-hand side (RHS) of eqn. (7) below. This would result in 
a steeper Emdt{Tc) (Jordan 2000), since it adds a positive 
term to the RHS of eqn. (10) below. An enthalpy flux is 
not included in either our full numerical solutions or the 
analytical approximations, because large enough systematic 
flows are not usually observed in the upper transition region 
and inner corona. Hence 

dFc(rc) ^ df,(re) 
dh ^ dh ' 

where the radiation losses are given by 

dF,(Te) _ 0.86Pi'P,ad(Te) 



(7) 



dh 



(8) 



and where Prad(T'c) is the radiative powerloss function. 
Or, using eqn. (3) for dh/dTe, 

dfr(ra) ^ V2gmdt(re)Prad(Te) 
dTe Te 

Eqn. (6) then becomes 

dlogSmdt(Te) _ ^ , 2dlogPo 2Emdt{T,fP,,i{n) 

- •^/^+ ,ji„„T^ „ „„ „._3/2 -y^^) 



(9) 



dlogTe 



dlogTe 



OMnPiTe^'^ 
1/2 



Here we approximate Prad(Te) by aTe" , where a is a con- 
stant (taken as 2.8 x 10"^® erg cm^ s~^ K^/^), on the 
grounds that we arc considering only coUisionally excited 
lines where 2 x 10^ K < To < lO'^ K. At higher temperatures 
continuum processes cause an increase in Prad(Te). Alter- 
natively, numerical values can be used, for example from 
CHIANTI (v6) (Dere et al. 2009). The form adopted here 
is useful in elucidating the physics since it results in the 
simplest analytical relations. 

The variation of the total pressure (including a non- 
thermal pressure term) with Te is included in our numerical 
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solutions using hydrostatic equilibrium. In the present pa- 
per, the non-thermal pressure is neglected, since in e Eri, 
above Te = 2 x 10 ' K, it does not exceed 0.08 of the gas 
pressure (Sim, 2002). Hydrostatic equilibrium should be a 
good approximation in the upper transition region and in- 
ner corona, since flows with velocities approaching the sound 
speed are not usually observed in these regions. 

The third term in eqn. (10) is given by hydrostatic equi- 
librium, and can be expressed as. 



.dlogPe 



, y/2nmiiEmdt{T,)g^T,Rl 



(11) 



dlogTe 0.86kBP|(ii. + /i)2 

where p is the mean molecular weight, taken to be 0.619, and 
Rt and are the stellar radius and surface gravity, respec- 
tively. The variation in this term is small in the mid tran- 
sition region, but becomes more important as Tg increases, 
since Eindt{Tc)Tc / increases with Te. 

Although wc drop the variation of Pc with Tc in general, 
we do use the integral of dPe/dTe to relate Pq and Pc, using 



2 2y2/^mHfl * f^" 
" 0.86kB 



To find an analytical solution (see Section 2.2) the varia- 
tion of the gravity with height has to be neglected. The full 
numerical solutions available for dwarf stars (e.g. Philippi- 
des 1996; Sim 2002; Ness & Jordan 2008) show that this 
variation is not very large in G/K dwarfs. E.g., in the final 
numerical model by Ness & Jordan (2008), using a radial 
height of < 3000 km at logTc = 5.30, gives g{ro)/g* > 0.99 
and g{rc)/g, > 0.90 at logTc = 6.53. 

From eqn. (10), ignoring the variation in pressure, a 
minimum in Erndt{Tc) occurs at some To, when 



Emt{To) 



3 X 0.86k 



4a 



(13) 



Thus at a given To, Emt{To) oc Po. 



2.2 Global constraints and resulting scalings 

At the top of the transition region/base of the corona, Pc(Tc) 
is the energy conducted back from the overlying heated 
corona. At To the conductive flux is Pc(To). The easiest way 
to use the global constraint that the overall net conductive 
flux is balanced by the total radiation losses is to use the 
approach of Rosner, Tucker & Vaiana (1978), but in the 
present work no explicit boundary conditions on the values 
of Pc(Tc) and Fc{To) have been imposed. 
From eqn. (7) one can write 



/: 



Pc(Te)dFc(Te) 



Pc(Te) 



dPr dh 
dh dTe 



dTe 



Hence, using eqns (4) and (8), 



Pc(To) 



: 0.86aK 



(14) 



(15) 



Then, on the left-hand side of eqn. (15), eqn. (3) can be used 
to express dTe/dh in terms of Emdt{Tc). On the RHS of eqn. 
(15) the pressure term is taken to be constant at Pq . This is 
close to the mean value of Pg in the full numerical solutions. 
The pressure variation between To and Tc is included in the 
conductive flux terms. The result is 



0.86k T^P^ 

A V 



:aPo'(Te-To) . 



(16) 



, ^£mt(Tc)2 Emt{To)^' 

Substituting for Emt{To) from eqn. (13), eqn. (16) can 
be rearranged to give 



Emt{Tc) 



0.86k TcP.^ 



4a Po (1 - 2^)1/2 ' 



(17) 



where Tc can be replaced by a general Te (> To), and sim- 
ilarly, Pc can be replaced by a general Pe, to give a general 
equation for EmdtiTc). 

Emt{Tc) can also be found by making the approxima- 
tion that the coronal emission is formed mainly over the 
first pressure-squared isothermal scale-height. This is jus- 
tified by using the equation of hydrostatic equilibrium to 
show that P^ decreases exponentially with height according 
to Pc exp[— (/i — hc)/H], where he is the height at the base 
of the corona and H = kBTc/2fj,m.ug*. Solar observations by 
Gibson et al. (1999) confirm the hydrostatic decrease in 
under the near-isothermal conditions in the low corona. In 
this case. 



Smt(Tc) 



0.86kB Pc 



(18) 



2/imH g*Tc' 

Combining eqns (17) and (18) allows Pc to be eliminated 
leading to an explicit expression for Po, 



0.86a ke (1-^)1/2 



(19) 



Thus Po scales as g^T^ , with only a weak dependence on 

To/Tc, and provided 5, is known, can be found from the 
values of To and Tc that match the observed behaviour of 
Pmdapp(Tc). 

Once Po is known, Emt{To) is known from eqn. (13) 
and can also be expressed in terms of g», Tc and To, i.e. 



£mt(To) = VS 



2akB (1 - |if )V2 
The ratio of Emt(Tc) to Emx(To) is given by 

gmt(Tc) ^ Tc Pc' 1 
Smt(To) To Pi 73(1 - |S.)V2 ■ 



(20) 



(21) 



To separate Emt{Tc) and Pc requires eqn. (12), which can 
be written as 



p2 

O 

P? 



1 + 



V2 J^^Emdt{Te) dTe 



(22) 



£;mt(Tc)Tc 

The detailed numerical models show that Po/Pc varies 
slowly when Tc is varied, because the last term in eqn. (22) 
is almost constant. The general form of EmdtiTc), when Tc 
in eqn. (17) is replaced by Te, can now be applied. This 
time, guided by the full calculations, we take the variable 
Pc' as Pc' and remove it from the integral. The integral of 
Te/[1 - (2/3)(To/Te)]^/' has a standard solution (see e.g. 
Jeffrey 1995), which reduces to 



8T^ 
9 



(1-F2) 



-dY 



To y(5 - 3y' 

9 (1-F2)2 



6 



■In 



\Y 



\Y + 1\ 



,(23) 



where F' = 1 - (2To/3Te). The first term on the RHS of 
eqn. (23), evaluated at Tc, dominates and reproduces the full 
solution to within 0.1 per cent. Thus, when the constants 
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in the general expressions for Emdt{Te) and Emt{Tc) are 
included, the last term in eqn. (22) becomes 



V2 J^^ Em(k{Te)dn 
Smt(Tc)Tc 

Hence Po/Pc can be found from cqns (22) and (24), and P^ 
can then be found, using Po from eqn. (19). 

Using eqn. (18) Emt{Tc) can then be expressed as 



Smt(Tc) 



2kBa (l-|^)[l + ^(i_|^)(l + |^)] 



.(25) 



This gives the scaling with g*T^ , expected from dimensional 
arguments, and also the absolute value, which arises from 
the energy balance adopted. The terms in To/Tc vary by 
only 1 per cent over the range of Tc in the models discussed 
below. Thus, for a given choice of k and a, only g^,^ Tc and 
To are required to evaluate Emt{Tc). 

We now make comparisons with results from Hcarn's 
(1975, 1977) formulation for minimum energy loss (mcl) 
coronae. Hearn applied dFc = —dFr to a corona, but took a 
partial differentiation with respect to the coronal tempera- 
ture, at constant coronal pressure. He assumed that To and 
Fc{To) can be neglected. 

Making the same assumptions, by dropping the second 
term in eqn. (15) and using eqn. (18), it can be seen that 



Po(mel) = 



0.86Q ke 



(26) 



However, in our approach, if Fc{To) tends to zero, then 
Emt{To) tends to infinity. Thus eqn. (16), together with eqn. 
(18) can be used to define a critical (ma^ximum) value of Po. 
This pressure cannot be reached when there is a minimum 
in Emdt{Te) at To. Po(crit) is given by 



Po(crit) 



K limn 



0.86a kB (1 - Si) 1/2 



(27) 



This approaches the value of Po(mel) when To is much 
smaller than Tc. If observations of density /pressure sensi- 
tive lines show clearly that Po exceeds Po(crit), then models 
including a fractional emitting area, or heating, additional 
to that provided by thermal conduction, must be considered. 

Although the difference between our predicted values of 
Po and Po(mel) arc not large, the advantage of our solutions 
is that they allow for the difference between Po and Pc in 
hydrostatic equilibrium and predict the values of Emt{To), 
as well as Emt{Tc). Hearn (1975, 1977) assumed a constant 
pressure with a value defined at the 'base of the corona'. 
Applying the pressure from eqn. (26) to find EnitiTc) from 
eqn. (18) gives larger values than found from our analytical 
solutions. 



3 RESULTS AND COMPARISONS WITH 
FULL MODELS 

In optimizing our full solution for the Emdt{Te) and 
Emda,pp{Te) of e Eri (Ness & Jordan 2008), we ran five mod- 
els with different values of Tc. The results for the full models 
are those for which a minimum in Emda,pp{Te) is just possi- 
ble at the chosen value of To = 2 x 10® K. Such solutions are 
found by gradually increasing the value of Eme,pp{Tc). As 



Table 1. Comparison of parameters from the numerical models 
for e Eri (lines labelled 'n') and from the analsrtical equations 
(lines labelled 'a'). The quantities involved, including their units 
are: log[Te(K)], log[Pe(cm-3 K)] and log[£;mt(re)(cm-5)]; logTo 
is fixed at 5.30. 



logTc 




6.50 


6.53 


6.55 


6.60 


6.65 


logPo 


(n) 


15.91 


15.97 


16.01 


16.10 


16.20 


logPo 


(a) 


15.86 


15.92 


15.96 


16.06 


16.16 


logEmt(ro) 


(n) 


27.50 


27.56 


27.60 


27.69 


27.79 


log£;mt(To) 


(a) 


27.37 


27.43 


27.47 


27.56 


27.66 


21og(Po/Pc) 


(n) 


0.254 


0.254 


0.253 


0.252 


0.251 


21og(Po/Pc) 


(a) 


0.236 


0.235 


0.235 


0.235 


0.235 


logPc 


(n) 


15.79 


15.84 


15.88 


15.98 


16.07 


logPc 


(a) 


15.75 


15.81 


15.85 


15.95 


16.04 


logEmt(rc) 


(n) 


28.22 


28.31 


28.36 


28.51 


28.66 


logEmt(rc) 


(a) 


28.10 


28.19 


28.25 


28.40 


28.55 



stressed earlier, the full numerical models arc in hydrostatic 
equilibrium, including a non-thermal pressure term based 
on observed line widths, and adopt a spherically symmetric 
atmosphere. Although we do not expect exact agreement 
between the analytical and the full models, it is of interest 
to examine the size of the differences between them, since 
the analytical models can be useful in the process of homing 
in on the optimum solution for a given star. 

Table 1 gives the values of the parameters discussed in 
the previous section, with those from the above full rmmeri- 
cal models given in the lines labelled 'n' and those from the 
analytical predictions in the lines labelled 'a'. The order of 
the parameters listed reflects the order in which the calcu- 
lations can be made, i.e., Po from eqn. (19), Emt(To) from 
eqn. (13), {Po/Pcf from eqns (22) and (24), Pc from eqns 
(19), (22) and (24), and Emt{Tc) from using Pc^ in eqns (17) 
or (18). 

The full numerical models give the same scaling laws 
for Po and Pc as found from the analytical approach, but 
with multiplication factors of about 1.1. Similarly, the full 
numerical emission measures [Emt{To) and Emt{Tc)] follow 
the same scaling laws as fotmd from the analytical approach, 
but are systematically larger by mean factors of about 1.3. 
The differences arise from the approximations to the elec- 
tron pressure used in the analytical equations, including the 
neglect of the non-thermal pressure term, and to the differ- 
ent geometries adopted. The variation of g with the radial 
distance r is also included in the rmmcrical solutions. 

Table 2 gives the combinations of parameters that ap- 
pear in eqn. (21) and in the fourth term of eqn. (10). 

It can be seen that the ratio of the emission measures at 
Tc and To given by eqn. (21) agrees better with the results 
from the full models than do the absolute values. The numer- 
ical models show that the ratio Emt{Tc) / PcTc is almost con- 
stant, as expected from the scaling laws, Emt{Tc) oc g{ro)To 
and Pc oc S'(rc)Tc^ 

In the rmmcrical solutions for e Eri, the origins of the 
values of the constants of proportionality in the empirical 
scaling laws are hard to pin down. The analytical expres- 
sions give similar scalings, but now the actual value of the 



6 C. Jordan, J.-U. Ness, S. A. Sim 



Table 2. Values of log[Emt{Tc)/Emt{To)] and Emt(Tc)/PcTc 
calculated from the full numerical models for e Eri (lines labelled 
'n') and the anals^ical solutions (lines labelled 'a'). Units as in 
Table 1. 



logTc 6.50 6.53 6.55 6.60 6.65 



log[£mt(Tc)/Bmt(To)] 


(n) 


0.72 


0.74 


0.76 


0.82 


0.87 


\og[Emt{Tc)/Emt(To)] 


(a) 


0.73 


0.76 


0.78 


0.83 


0.88 


log[Emt(rc)/PcTc] 


(n) 


5.93 


5.93 


5.93 


5.93 


5.94 


iog[£mt(rc)/Pcrc] 


(a) 


5.85 


5.85 


5.85 


5.85 


5.85 



Table 3. Predicted analytical values of Po, Pc, Emt(To) and 
Emt{Tc) for Procyon, together with combined parameters. Re- 
sults using log To = 5.5 are given in the upper part of the table, 
and those using log To = 5.7 in the lower part. Units as in Table 
1. 



logTc 


6.20 


6.25 


6.30 


log To = 5.5 








logPo 


14.62 


14.71 


14.81 


logPc 


14.50 


14.59 


14.69 


log£;mt(To) 


26.32 


26.41 


26.51 


logEmt(Tc) 


26.57 


26.71 


26.86 


log[Bmt(Tc)/PcTc] 


5.87 


5.87 


5.87 


log To = 5.7 








logPo 


14.64 


14.73 


14.82 


logPc 


14.52 


14.61 


14.70 


\ogEmt(To) 


26.54 


26.63 


26.73 


logSmt(Tc) 


26.61 


26.75 


26.89 


log[Emt(Tc)/PcTc] 


5.89 


5.89 


5.89 



et al. (2002) also analyzed spectra from the LETGS and 

the XMM-Newton satellite, but made a global 3-T fit us- 
ing the SPEX computer package, rather than an individual 
line-based approach. Wood et aJ. (1996) analyzed the line 
widths and redshifts measured from spectra obtained with 
the GHRS instrument on the HST. The line widths were 
found to bo larger than in the main-sequence stars, such 
as e Eri (Sim & Jordan 2003b). In summary, in Procyon, 
log[To(K)] lies between 5.5 and 5.7 and log[rc(K)] lies be- 
tween 6.2 and 6.3. The equation of hydrostatic equilibrium 
should strictly include the effects of any non-thermal pres- 
sure associated with the larger line widths observed. 

We adopt the following stellar properties: a distance 
of 3.53 pc (Girard et al. 2000), an angular diameter of 
5.51 mas (Mozurkewich et al. 1991) and hence R* = 2.09Rq ; 
a mass of M» = I.5M0 (Girard et al. 2000) and hence 
log[g.(cm s"^)] = 3.98. 

The results of applying the analytical expressions given 
in Section 2 are summarized in Table 3, for log[To(K)] = 5.5 
and 5.7, and log[Tc(K)] = 6.2, 6.25 and 6.3. 

As yet there are no completely satisfactory numeri- 
cal models of the outer atmosphere of Procyon. Philppides 
(1996) used observations from ROSAT {ROentgen SATellit) 
to find Emapp{Tc) and Tc, and made models in a spheri- 
cally symmetric geometry, including a non-thermal pressure 
term. She noted that the latter term causes Pc to increase 
with Tc within the transition region, before decreasing again 
by Tc. The later observations with EUVE (Drake et al. 1995; 
Schmitt et al. 1996) and both the LETGS and EUVE (Sanz- 
Forcada et al. 2004) showed that, as expected, the earlier 
1-temperature and 2-tempcrature fits to ROSAT spectra 
overestimated the coronal emission measure. Sim (1998, un- 
published MPhys project; 2002) made a model of the chro- 
mosphere and lower transition region using line fluxes and 
widths measured by Wood et al. (1996), which supersede the 
fluxes from Copernicus by Jordan et al. (1986), and included 
the radiative transfer in lines formed up to ~ 2 x 10* K. He 
adopted a plane parallel atmosphere up to 3 x 10^ K and a 



constant of proportionality can be clearly tracked back to 
the global energy balance equation assumed. 



4 EXAMPLE APPLICATION 

Here we apply the analytical expressions to Procyon (HD 
61421, F5 IV-V) to illustrate what can be learnt before de- 
tailed modelling is carried out. Earlier work on Procyon has 
shown that it is difficult to reconcile different measurements 
of Pc without invoking limited areas of emission or the pres- 
ence of active region loops (Schmitt et al. 1985, 1996; Jor- 
dan et al. 1986). Because of the lower gravity and larger line 
widths in Procyon, compared with those of cool dwaxf stars, 
the simple methods used in Section 2 are expected to be less 
accurate than for e Eri. For a quiescent corona, the nature of 
the heating flux is not relevant in our energy balance model, 
so either an MHD wave flux or an acoustic wave flux, as 
suggested by Mullen & Cheng (1994), is possible. 

Emission line fluxes are available for Procyon from a 
number of the instruments mentioned in Section 1. Jordan 
et al. (1986) analyzed spectra obtained with the lUE and the 
Einstein Observatory and derived emission measures, but 
not a mean emission measure distribution. Early observer 
tions with Copernicus were used to constrain the line emis- 
sion measure around 3 x 10 ' K. They found that the ratio of 
the coronal emission measure to that of the lower transition 
region was significantly smaller than in main-sequence dwarf 
stars and that the mean coronal temperature was around 
1.5 X lO** K. This was in broad agreement with earlier work 
by Schmitt et al. (1985), who used data from the Einstein 
Observatory. In particular, even with the uncertainty in the 
line fluxes from Copernicus, it appeared that To > 3 x 10® K, 
rather than 2 x 10^ K in the main-sequence stars. Drake 
et al. (1995) included data obtained with the EUVE, plus 
some adjusted data from Copernicus, to produce the emis- 
sion measure distribution above Tg ~ 1.6 x 10® K, while 
Sanz-Forcada, Brickhouse & Dupree (2003) combined the 
observations from the lUE and the EUVE to improve the 
overall Emda,pp{Tc). Although there were differences in de- 
tail, owing to the abundances and atomic data adopted, the 
Smdapp(Te) found by Sanz-Forcada et al. (2003) showed 
a similar form to that indicated in Jordan et al. (1986), 
but with To ~ 4 X Uf K. Sanz-Forcada et al. (2004), im- 
proved the higher temperature part of the EMD using spec- 
tra obtained with the LETGS on the Chandra satellite. 
They found To ~ 5 x 10® K and Tc ~ 2 x 10® K. Raassen 
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Table 4. Values of Po, Pc, Smapp(To) and -Emapp(Tc), from 
numerical models and observations. Units as in Table 1. 



log Tc 6.30 



log To = 5.5 

log Po 14.81" 
logPc 14.73" 
logEmapp(To) 26.27'' 
logEmapp(rc) 26.84^= 
27.16'' 



" Prom numerical model by Sim (2002). 

Prom observations fitted by Sim (1998). 
° Starting value adopted by Sim (1998). 

Prom observations fitted by Sanz-Porcada et al. (2004), scaled 
as described in the text. 

spherically symmetric atmosphere at higher temperatures, 
including the non-thermal pressure term throughout the at- 
mosphere. However, he interpolated Emda,pp{Te) between 
~ 3 X 10'^ K and a coronal temperature of 2 x 10® K, rather 
than making an energy balance model. 

Here we use the work by Sanz-Forcada et al. (2004) to 
give information on Emdi,pp{Te), but adopt an iron abun- 
dance of log A^Fc = 7.51, rather than their value of 7.32 (on 
the scale where log A^h is 12.0). Their absolute scale has been 
changed to match our definition of Errit (0.3) given by eqn. 
(1). Their maximum value of the apparent volume emission 
measure at log[Tc(K)] = 6.30 has been divided by to 
find i?mapp(T'c). Because of this conversion, in a spherically 
symmetric geometry the apparent and true emission mea- 
sures are related by 

Em^pp(n) = f{r)G(r)Emt(n) , (28) 

where /(r) = {r/R, f and G(r) = 0.5(1 + ^[1 - (l//(r))]) 
is the fraction of emitted photons not intercepted by the 
star. Thus in an extended atmosphere the product f(r)G(r) 
can be greater than 1, so that iJmapp(rc) can be larger than 
Smt(Tc). 

Table 4 gives the observed values of i5mapp(T'o) and 
Smapp(Tc) from Sim (1998) and Sanz-Forcada et al. (2004), 
respectively. Values of log Po and log Pc from Sim's (2002) 
numerical models are also given. 

We can now compare the predicted values of EnitiTo) 
and Emt{Tc), given in Table 3 for log[Tc(K)] = 6.30, with 
the values of Smapp(To) and Ema,pp{Tc) given in Table 4. 
The value of log £;mapp(To), when log[To(K)] = 5.50, found 
by Sim (1998) is smaller than the value of logi5mt(ro). In 
a plane parallel atmosphere these are expected to differ by 
about a factor of two, which is indeed the case. Bearing in 
mind the comparisons between the analytical and numeri- 
cal results for e Eri, this suggests that the emission at this 
height is not from a highly restricted area of the atmosphere. 
The value of Emt{Tc) is similar to, or less than the value 
of _Emapp(Tc). Using the pressure-squared isothermal scale 
height (H) in the corona, plus P* , to estimate r at Tc leads 
to /(r)G(r) ~ 0.85. Given the uncertainties in, and diS'er- 
ences between the atomic data used by the above authors, 
the disagreement between the predicted and observed values 
is not large and is in the direction found from the numerical 



and analytical results for e Eri. However, Procyon is a case 
where full numerical models are required to make more de- 
tailed comparisons between observed and predicted results. 

It is very difficult to establish values of A'e in Pro- 
cyon. Schmitt et al. (1996) made a careful analysis of den- 
sity sensitive lines of Fex to xiv observed with the EUVE 
and concluded that the average value of A'o at coronal tem- 
peratures lies between 10^ and 10^° cm~^, with a value of 
3 X 10^ cm~^ being adopted. However, for some ions the re- 
sults depended on weak lines, and a range of A^e was found 
from different pairs of lines within a given ion. As Schmitt et 
al. (1996) point out, the uncertainties in iVe mask any sys- 
tematic variation of A^e or Pc with the stage of ionization. 
At around log[Tc(K)] = 6.20 — 6.30, their smallest value of 
log[Pc(cm^^ K)] = 15.2 is significantly larger than those 
given in Table 3. 

Ness et al. (2001) measured the ratio of the fluxes in 
the forbidden line (Is^ ""^So - ls2s '^Si) and intcrsystcm (plus 
quadrupolc) line (Is'^ ^So - ls2p '^Pi,2) in the Hcl-like ions 
C V, N VI and O VII, using spectra obtained with the LETGS. 
We have found revised flux ratios by also including further 
spectra available from more recent observations with the 
LETGS. We have also updated the values of (and limits 
on) A'"e by using CHIANTI (v6) (Dere et al. 2009). In the 
Sun, only the flux ratio in Cv is afiected by the photo- 
spheric/chromospheric radiation field at the wavelength of 
the *Si - ^P transitions (Gabriel & Jordan 1969). In Pro- 
cyon, because of the higher pliotosphcric temperature, Ness 
et al. (2001) found that this photoexcitation is also signifi- 
cant in Nvi. As a result, A'e is not constrained by the ob- 
served flux ratio in Cv. In Nvi, the combined uncertainties 
in the radiation field and the measured flux ratio are too 
large to yield a definitive value of log[Afc(cm^'')], which can 
lie between < 8.0 and 10.0. Ovii provides the best diag- 
nostic, since the effects of the radiation field are small and 
make little difference to the value of log[A^c(cm~^)] = 9.13 
obtained (when these effects are included). But the uncer- 
tainty of about ± 10 per cent in the measured flux ratio 
of 3.60 includes the value of 3.75 at log[Afe(cm-^)] = 8.0. 
The observed ratio yields log[Pe(cm~^ K)] = 15.43, with an 
upper limit of 15.99 and a lower limit of < 14.3. 

Liang, Zhao & Shi (2006) investigated density-sensitive 
X-ray lines of Si X, using LETGS fluxes from Raassen et al. 
(2002). From the strongest lines at 50.524 and 50.691 A they 
found log[Are(cm"^)] ~ 8.41-8.45 at Te = 1.26x 10® K, and 
that calculations by other authors gave only slightly higher 
values. The resulting values of log[Pe(cm~^K)] were in the 
range 14.5 - 14.7. These are similar to the value found from 
Sim's (2002) numerical model and to the value predicted by 
the analytical scaling laws. The Six lines are expected to be 
formed at the same Tc as those of N VI. However, in Si x, the 
high density limit for the flux ratio occurs at a lower value 
of A'e than in Nvi. If the emitting regions contained both a 
quiescent corona and active regions (of smaller area), then 
the higher density regions could contribute relatively less to 
the electron density measured from Six. Thus one cannot 
rule out the possibility that N vi could detect a higher value 

of ATe. 

Jordan et al. (1986) used simple line opacity arguments 
to deduce that at Tc = 2 x 10^ K, log[Pc(cm-^ K)]< 14.4. 
This is a factor of two lower than the value of 14.76 in the 
model by Sim (2002). Wood et al. (1996) made a more so- 
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phisticated estimate of line opacity efTects and deduced that 
the profile of the line of Sim at 1206 A implied a pressure of 
log[Pc(cm~^ K)]< 14.8 at log[Tc(K)] = 4.70, in reasonable 
agreement with the value of 14.68 in Sim's (2002) model. 
However, the intersystem lines of O iv analyzed by Wood et 
al. (1996) led to log[Pe(cm-3 K)]~ 15.0 (or < 15.5, when 
possible errors were considered). But as they pointed out, 
the O IV flux ratios were not far from those expected in the 
low-density limit. Analyses of spectra observed recently by 
Ayre^ with the STIS should give at least improved limits 
on values of Pc- The values of log Po from the methods used 
here are given in Table 3 and agree well with Sim's (2002) 
model value of 14.81. 

Overall, the higher coronal values of Pc found by 
Schmitt et al. (1996), and from our present analysis of Ovil, 
could be reconciled with the values from the analytical pre- 
dictions and the numerical models if higher pressure active 
region loop structures were present, as well as a quiescent 
corona (see also Schmitt et al. 1985, Jordan et al. 1986). 
Such active regions could also contribute to Emd^pp near 
its apparent maximum value. However, it seems very un- 
likely that the quiescent atmosphere has pressures that sig- 
nificantly exceed those given in Tables 3 and 4. Whether or 
not the upper part of the quiescent atmosphere is heated 
by acoustic or MHD waves is still an open question, but 
in the time-averaged acoustic heating model by Mullan & 
Cheng (1994), the maximum coronal temperature is only 
6.5 xlO® K. If this temperature is adopted for Tc, the pre- 
dicted value of Emt{To) is about a factor of 5 smaller than 
the observed value given in Table 4, so any acoustic heating 
must lead either to a hotter corona, or be limited to the 
region below To. 

By applying the analytical solutions we have found that 
the value of Emt{To) is acceptable without invoking a lim- 
ited area of emission and that current acoustic heating mod- 
els are not entirely satisfactory. We have also predicted the 
values of Po in the quiescent transition region and corona. 
While the former is acceptable, higher pressures cannot be 
ruled out above about Te > 10^ K. Procyon is clearly a star 
for which full numerical models are required. These will be 
carried out following the methods used by Ness & Jordan 
(2008) in studies of e Eri, and will include recent atomic 
data for the emission lines used. 



5 DISCUSSION AND CONCLUSIONS 

Full numerical solutions giving the EmdtiTc) and 
Emd^pp(Te) of e Eri were carried out by Ness & Jordan 
(2008). The solutions for different values of Tc and fixed To 
and gifc) showed scaling laws between these parameters and 
Po, Pc, Emt{To) and Emt{Tc). These solutions were calcu- 
lated in spherical symmetry and in hydrostatic equilibrium, 
including a non-thermal pressure term. The link between 
the scaling laws and the assumed energy balance was not 
obvious. Here, analytical solutions for a plane parallel at- 
mosphere, with some approximations for the variation of Pc 
are presented, which give scaling laws that can be linked di- 
rectly to the energy balance adopted. These analytical solu- 
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tions reproduce well the values of Po and Pc and give values 
of ErritiTo) and Emt{Tc) that are smaller than in the full 
solutions by less than about a factor of 1.3. The analytical 
solutions are therefore useful in finding the best initial con- 
ditions in full numerical solutions. Only the values of g*. To 
and Tc are required, the latter two being simple to measure. 

When a 'critical solution' is found by setting the con- 
ductive flux at To to zero, and when To is very much less 
than Tc, our value of Po(crit) is the same as that given by 
Hearn's (1975, 1977) minimum energy loss hypothesis. In 
our approach, the use of the condition that Emdi{TB) passes 
through a minimum at some observed To has the advan- 
tage that Pc and Emt(To) can be found, as well as Po and 
-Bmt(rc). 

We have applied our results to Procyon, for which the 
approximations made should be less accurate than for e Eri. 
It is difficult to measure the electron density in Procyon, 
but our results rule out suggestions that in the spatially 
averaged atmosphere the values significantly exceed about 
7 X 10^ cm~"^ in the transition region above 10^ K, or about 
3 X 10* cm~^ in the quiescent corona. Improvements to the 
density-sensitive line ratios in the EUVE wavelength range 
will have to await the flight of new instruments. 

Some earlier problems in reconciling coronal and tran- 
sition region pressures have been reduced by using the full 
Emda.pp{Tc) , rather than single Tc fits to earlier measure- 
ments of X-ray fluxes. If the higher electron densities found 
from the EUVE and possibly from the LETGS could be 
confirmed, these would show that active region material is 
indeed present (Schmitt et al. 1985; Jordan et al. 1986). 
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